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ABSTRACT 

Context. Debris discs around main-sequence stars indicate the presence of larger rocky bodies. The components of the nearby, solar-type binary 
a Centauri have metallicities that are higher than solar, which is thought to promote giant planet formation. 

Aims. We aim to determine the level of emission from debris around the stars in the a Cen system. This requires knowledge of their photospheres. 
Having already detected the temperature minimum, T min , of a Cen A at far-infrared wavelengths, we here attempt to do the same for the more 
active companion crCenB. Using the a Cen stars as templates, we study the possible effects that T m \ n may have on the detectability of unresolved 
dust discs around other stars. 

Methods. We used Herschel- PACS, Herschel-SPIRE, and APEX-LABOCA photometry to determine the stellar spectral energy distributions in the 
far infrared and submillimetre. In addition, we used APEX-SHeFI observations for spectral line mapping to study the complex background around 
a Cen seen in the photometric images. Models of stellar atmospheres and of particulate discs, based on particle simulations and in conjunction 
with radiative transfer calculations, were used to estimate the amount of debris around these stars. 

Results. For solar-type stars more distant than a Cen, a fractional dust luminosity / d = E dust /L star ~ 2 x 10“ 7 could account for SEDs that do 
not exhibit the T min effect. This is comparable to estimates of / d for the Edgeworth-Kuiper belt of the solar system. In contrast to the far infrared, 
slight excesses at the 2.5cr level are observed at 24 /um for both a Cen A and B, which, if interpreted as due to zodiacal-type dust emission, would 
correspond to / d ~ (1-3) x 10~ 5 , i.e. some 10 2 times that of the local zodiacal cloud. Assuming simple power-law size distributions of the dust 
grains, dynamical disc modelling leads to rough mass estimates of the putative Zodi belts around the a Cen stars, viz. ^4x 10 -6 M ^ of 4 to 1000 /um 
size grains, distributed according to n(a ) oc a~ 35 . Similarly, for filled-in T m \ n emission, corresponding Edgeworth-Kuiper belts could account for 
~10 _3 M(£ of dust. 

Conclusions. Our far-infrared observations lead to estimates of upper limits to the amount of circumstellar dust around the stars a Cen A and B. 
Light scattered and/or thermally emitted by exo-Zodi discs will have profound implications for future spectroscopic missions designed to search 
for biomarkers in the atmospheres of Earth-like planets. The far-infrared spectral energy distribution of a Cen B is marginally consistent with the 
presence of a minimum temperature region in the upper atmosphere of the star. We also show that an a Cen A-like temperature minimum may 
result in an erroneous apprehension about the presence of dust around other, more distant stars. 

Key words, stars: individual: Alpha Centauri - binaries: general - circumstellar matter - infrared: stars - infrared: planetary systems - 
submillimeter: stars 


1. Introduction 

The a Centauri system lies at a distance of only 1.3 pc ( n - 
747 .1 ± 1.2 mas, Soderhjelm 1999), with the G2 V star a Cen A 
(HIP 71683, HD 128620) often considered a solar twin. Together 
with the K 1 star crCenB (HIP 71681, HD 128621), these stars 
are gravitationally bound in a binary system, with an orbital pe- 
riod of close to 80 years and a semi-major axis (24 AU) midway 
between those of the planets Uranus and Neptune in the solar 
system. A third star, Proxima Centauri, about 2° southwest of the 
binary shares a similar proper motion with them and seems cur- 
rently to be bound to trCenAB, although the M6 star trCenC 
(HIP 70890) is separated by about 15 000 AU. 


* Based on observations with Herschel which is an ESA space ob- 
servatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 

** And also based on observations with APEX, which is a 12 m di- 

ameter submillimetre telescope at 5100 m altitude on Llano Chajnantor 
in Chile. The telescope is operated by Onsala Space Observatory, Max- 
Planck-Institut fur Radioastronomie (MPIfR), and European Southern 
Observatory (ESO). 


The question as to whether there are also planets around our 
solar-like neighbours has intrigued laymen and scientists alike. 
The observed higher metallicities in the atmospheres of a Cen A 
and B could argue in favour of the existence of planets around 
these stars (Maldonado et al. 2012, and references therein). The 
proximity of aCe n should allow for highly sensitive observa- 
tions at high angular resolution with a variety of techniques. 

We know today that binarity is not an intrinsic obstacle to 
planet formation, since more than 12% of all known exoplan- 
ets are seen to be associated with multiple systems (Roell et al. 
2012). Even if most of these systems have very wide separations 
(>100 AU) for which binarity might have a limited effect in the 
vicinity of each star, a handful of planets have been detected in 
tight binaries of separation ~20 AU (e.g., y Cep, HD 196885), 
comparable to that of crCen (Desidera & Barbieri 2007; Roell 
et al. 2012). The presence of these planets poses a great chal- 
lenge to the classical core- accretion scenario, which encounters 
great difficulties in such highly perturbed environments (see re- 
view in Thebault 2011). 

For the specific case of a Centauri, radial velocity (RV) ob- 
servations indicate that no planets of mass >2.5 Mj U pi ter exist 
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inside 4AU of each star (Endl et al. 2001). For their part, 
theoretical models seem to indicate that in situ planet forma- 
tion is indeed difficult in vast regions around each star, be- 
cause the outer limit for planet accretion around either star is 
-0.5-0. 75 AU in the most pessimistic studies (Thebault et al. 
2009) and -1-1.5 AU in the most optimistic ones (e.g. Xie 
et al. 2010; Paardekooper & Leinhardt 2010). However, these 
estimates open up the possibility that planet formation should be 
possible in the habitable zone (HZ) of aCenB, which extends 
between 0.5 and 0.9 AU from the star (Guedes et al. 2008). Very 
recently, based on a substantial body of RV data, Dumusque 
et al. (2012) have proposed that an Earth-mass planet orbits 
or Cen B with a three-day period (but see Hatzes 2013). In other 
words, the radial distance of crCenBb (0.04 AU), which corre- 
sponds to only nine stellar radii, is evidently far inside the HZ, 
and therefore the surface conditions should be far from able to 
support any form of life as we know it. 

Based on sensitive Herschel (Pilbratt et al. 2010) observa- 
tions, a relatively large fraction of stars with known planets ex- 
hibit detectable far-infrared (FIR) excess emission due to cool 
circumstellar dust (Eiroa et al. 2011; Marshall et al. 2014; Krivov 
et al. 2013), akin to the debris found in the asteroid (2-3 AU) 
and Edgeworth- Kuiper (30-55 AU) belts of the solar system. As 
part of the Herschel open time key programme DUst around 
NEarby Stars (DUNES; Eiroa et al. 2013), we observed crCen 
to search for dust emission associated with the stars, which is 
thought to develop on planetesimal size scales and to be ground 
down to detectable grain sizes by mutual collisions. Persistent 
debris around the stars would be in discs of a few AU in size and 
would re-emit intercepted starlight in the near- to mid-infrared. 
To within 16%, ISO-SWS observations did not detect any ex- 
cess above the photosphere of crCenA between 2.4 and 12 /mi 
(Decin et al. 2003, and Fig. 7 below). On the other hand, a 
circumbinary Edgeworth- Kuiper belt analogue would be much 
larger and the dust much cooler, so that it would emit predom- 
inantly at FIR and submillimetre (submm) wavelengths. Such a 
belt would be spatially resolved with the Herschel beam, in prin- 
ciple allowing the detection of structures due to dynamical inter- 
actions with, say, a binary companion and/or giant planets (e.g. 
Wyatt 2008, and references therein), but its surface brightness 
could be expected to be very low, rendering such observations 
very difficult. 

The DUNES programme focusses on nearby solar- type stars, 
and the observations with Herschel - PACS (Poglitsch et al. 2010) 
at 100/rni and 160 /mi are aimed at detecting the stellar photo- 
spheres at an S /N > 5 at 100 /mi and has observed 133 nearby 
Sun-like stars (FGK type, d < 25 pc). Prior to Herschel , little 
data at long wavelengths that have high photometric quality are 
available for a Cen A and B. One reason for this is probably de- 
tector saturation issues due to their brightness (e.g. WISE bands 
W 1-W4), and another is due to contamination in the large beams 
by confusing emission near the galactic plane (e.g. IRAS, ISO- 
PHOT, and AKARI data). 

Due to their proximity and with an age comparable to that of 
the Sun (4.85 Gyr, Thevenin et al. 2002), a Cen is an excellent 
astrophysical laboratory for “normal” low-mass stars, otherwise 
known to be very difficult to calibrate, not the least with respect 
to their ages. From numerous literature sources, Torres et al. 
(2010) have compiled the currently best available basic stellar 
parameters of the a Cen system. The given errors on the phys- 
ical quantities are generally small. However, whereas the tab- 
ulated uncertainty of the effective temperature of e.g. crCenA 
is less than half a percent, the observed spread in Table 1 of 
Porto de Mello et al. (2008) corresponds to more than ten times 


this much. On the other hand, the radii given by Torres et al. 
(2010) are those directly measured by Kervella et al. (2003) us- 
ing interferometry, with errors of 0.2% and 0.5% for A and B, 
respectively (Bigot et al. 2006). Masses have been obtained from 
astro seismology and are good to within 0.6% for both compo- 
nents (Thevenin et al. 2002). 

For such an impressive record of accuracy for the stellar pa- 
rameters of the a Cen components, it should be possible to con- 
struct theoretical model photospheres with which observations 
can be directly compared to a high level of precision. Here, we 
report PACS observations of a Cen at 100 yum and 160 yum. These 
single-epoch data are complemented by LAB OCA (Siringo et al. 
2009) data at 870/rni obtained during two different epochs. 
The LAB OCA observations primarily address two issues: the 
large proper motion (377 yr -1 ) should enable the discrimination 
against background confusion and, together with SPIRE pho- 
tometry (see below), these submm data should also provide valu- 
able constraints on the spectral energy distributions (SEDs). This 
could potentially be useful for quantifying some of the proper- 
ties of the emitting dust and/or to gauge the temperature minima 
at the base of the stellar chromospheres. A clear understanding 
of the latter is crucial when attempting to determine extremely 
low levels of cool circumstellar dust emission. 

The paper is organised as follows. Section 2 outlines our 
observations with various facilities, both from space and the 
ground. Since we are aiming at low-level detections, the reduc- 
tion of these data is described in detail. Our primary results are 
communicated in Sect. 3. In the discussion Sect. 4, we examine 
the lower chromospheres of the a Cen stars in terms of the radi- 
ation temperatures from their FIR photospheres. Possible contri- 
butions to the FIR/submm SEDs by dust are also addressed, us- 
ing both analytical estimations and detailed numerical models. 
Finally, Sect. 5 provides a quick overview of our conclusions. 

2. Observations and data reduction 

2.1. Herschel 

In the framework of our observing programme, i.e. the DUNES 
open time key programme, PACS photometric images were ob- 
tained at 100 /urn and 160 /mi. In addition, from the Hi-GAL sur- 
vey (PI S. Molinari), we acquired archive data for a Cen at 70 jm i 
and 160 /mi obtained with PACS and at 250 /mi, 350 yum, and 
500 yum with SPIRE (Griffin et al. 2010). The relative position of 
the stars in the sky during the observational period is shown in 
Fig. 1. 

2.1.1. DUNES: PACS lOOyum and 160 yum 

PACS scan maps of crCen were obtained at 100 yum and 160 /mi 
at two array orientations (70° and 110°) to suppress detector 
striping. The selected scan speed was the intermediate setting, 
i.e. 20" s“\ determining the PSF at the two wavelengths (777 
and 12", respectively). The 100/rni filter spans the region 85- 
130 /mi and the 160 yum filter 130-2 10 /mi and the observations 
at 100 and 160/rni were made simultaneously. The data were 
reduced with HIPE v.8.0.1. The native pixel sizes are 372 and 
6'/4 at 100 and 160/rni, respectively, and in the reduced images 
the resampling resulted in square pixels of 1" at 100 /mi and 2" 
at 160 /mi. 

The two stellar components are approximated by model in- 
strument PSFs of the appropriate wavelength based on an ob- 
servation of a Boo rotated to match the telescope position an- 
gle of the a Cen observations. It is important to match the 
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Table 1 . Observing log. 


Instrument/ 

mode 

Obs/Pgm ID 

Wavelength 
T eff (yum) 

Beam width 
HPBW (") 

Observing date 
year-mo-day 

hnt 

(s) 

Centre coordinates" 
h m s ° ' " 

Offset" 

(") 

Herschel - PACS 

1342224848 

100 

777 

2011-07-29 

542 

14 39 30.200 -60 49 59.66 

1.7 


1342224849 

160 

1173 

2011-07-29 

542 

14 39 30.115 -60 49 59.54 

1.3 

SpirePacsParallel 

1342203280 

70 

579 x 1272 

2010-08-21 

9490 

14 28 59.831 -60 39 01.60 

2°x 2° map 


1342203281 

160 

1 1 76 x 1576 

2010-08-21 






250 

17'/6 C 

2010-08-21 






350 

23 "9 C 

2010-08-21 






500 

35'/2 c 

2010-08-21 




APEX-LABOCA 

384.C-1025(A) 

870 

1975 

2009-09-19 

7147 

14 39 31.725 -60 49 59.10 

3.7 


380.C-3044(A) 

870 

1 975 

2007-1 1-10 to 13 

34026 

14 39 32.349 -60 50 00.00 

3.6 

APEX-SHeFI 

090.F-9322(A) 

1300 

27'/ 1 

2012-08-16 

2646 

14 39 35.060 -60 50 15.1 

sp. line map 


Notes. (a) Observed equatorial coordinates (J2000) towards the target, i.e. the centre of the fitted point source, which refers to the primary a Cen A, 
except for Hi-GAL and APEX-SHeFI maps. Offset of observed with respect to intended (= commanded) position. Except for SHeFI, coordinates 

are corrected for the stellar proper and orbital motions. Not applicable (na) to survey data. (c) According to the SPIRE manual, see http : // 
herschel.esac.esa.int/Docs/SPIRE/html/, the average beam area is 423 arcsec 2 at 250yum, 751arcsec 2 at 350yum and 1587arcsec 2 at 
500 /mi. 



Fig. 1 . Binary a Cen AB as it appears in the sky on various occasions. 
Shown is the orbit of a Cen B with respect to the primary A, which is 
at the origin, and with north up and east to the left. During the observ- 
ing period, the separation of the stars was greatest in 2005 and smallest 
in 2011. The relative positions in the PACS (blue) and LABOCA (red) 
maps are indicated. The proper motion of a Cen for the time span be- 
tween the LABOCA and PACS observations is shown by the dashed 
arrow. The separation between the stars at the time of the MIPS obser- 
vations is indicated by the green dot. Orbital elements are adopted from 
Pourbaix et al. (2002). 

orientation of the PSF due to the non-circular tri-lobal structure 
of the Herschel - PACS PSF, which exists below the 10% peak 
flux level. 

Fitting of the two components was carried out by subtracting 
the components in series, starting with the brighter one. Each 
component model PSF was scaled to the estimated peak bright- 
ness and shifted to the required positional offset from the ob- 
served source peak before subtraction. The position offset of B 
relative to A was fixed. 


As described in detail by Eiroa et al. (2013), the level of the 
background and the rms sky noise were estimated by calculat- 
ing the mean and standard deviation of 25 boxes sized 9" x 9" 
at 100 //m and 14" x 14" at 160 yum scattered randomly at posi- 
tions lying between 30"-60" from the image centre and within 
the area of which no pixel was brighter than twice the standard 
deviation of all non-zero pixels in the image (the threshold crite- 
rion for source/non- source determination while high pass filter- 
ing during map creation). The calibration uncertainty was as- 
sumed to be 5% for both 100 //m and 160 /rm (Balog et al. 
2013) 1 . 

Aperture (and potential colour) corrections of the stellar flux 
densities and sky noise corrections for correlated noise in the 
super-sampled images are described in the technical notes PICC- 
ME-TN-037 and PICC-ME-TN-038 and the NHSC/PACS home 
page 2 . 

2.1.2. Hi-GAL: PACS 70yum and 160yum 

As part of the Hi-GAL program, these data were obtained at 
a different scan speeds, i.e. the fast mode at 60" s -1 . The 
field 314_0, containing crCentauri, was observed at both wave- 
lengths simultaneously and in parallel with the SPIRE instru- 
ment (see next section). The scanned area subtends 2° x2° and 
these archive data are reduced to Level 2.5. Compared to the 
data at longer wavelengths, background problems are much less 
severe at 70 /rm. The PACS 70 yum data provide the highest angu- 
lar resolution of all data presented here (Table 1). As before, no 
further colour correction was required, but the proper aperture 
correction (1.22) was applied. After re-binning, the pixel size is 
the same as that at 100/mi, viz. one square arcsecond, whereas 
the 160/rni pixels are, as before, two arcseconds squared. 

2.1.3. Hi-GAL: SPIRE 250yum, 350yum, and 500 yum 

At the relatively low resolution of the SPIRE observations, 
the strong and varying background close to crCen presented a 
considerable challenge for the flux measurements. However, for 

1 The PACS calibration scheme is described further in detail at 
http : //herschel . esac . esa . int/twiki/bin/view/Public/ 
PacsCalibrationWeb#PACS_instrument_and_calibration 

2 https : //nhscsci . ipac . caltech . edu/sc/index . php/Pacs/ 
AbsoluteCalibration 
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Table 2. Photometry and FIR/flux densities of a Centauri. 


4 e ff 

(jjm) 

a Cen A 
(Jy) 

a Cen B 
s v (Jy) 

Photometry 
& reference 

0.440 

2215 + 41 

536+ 10 

5(1) 

0.550 

3640 ± 67 

1050 + 19 

V(l) 

0.790 

4814 + 89 

1654 + 30 

/(l) 

0.440 

2356 + 43 

572+ 10 

6(2) 

0.550 

3606 + 66 

1059 + 20 

V (2) 

0.640 

4259 + 78 

1387 + 26 

R c(2) 

0.790 

4784 + 88 

1666 + 31 

I c(2) 

1.215 

4658 + 86 

1645 + 30 

7(3) 

1.654 

3744 + 69 

1649 + 31 

H{ 3) 

2.179 

2561 +47 

1139 + 21 

K( 3) 

3.547 

1194 + 22 

521 + 10 

L( 3) 

4.769 

592+ 11 

258 + 5 

M(3) 

24 

30.84 + 0.76 

13.63 + 0.33 

MIPS (4) 

70 

3.35 + 0.28 

1.49 + 0.28 

PACS (5) 

100 

1.41+0.05 

0.67 + 0.037 

PACS (6) 

160 

0.56 + 0.06 

0.21 + 0.06 

* PACS (5), (6) 

250 

0.24 + 0.05 

0.11+0.05 

* SPIRE (5) 

350 

0.145 + 0.028 

0.064 + 0.028 

* SPIRE (5) 

500 

0.08 + 0.03 

0.04 + 0.03 

* SPIRE (5) 

870 

0.028 + 0.007 

0.012 + 0.007 

* LABOCA (7) 


Notes. ( * } Asterisks indicate values determined using S v>a /Sv,b = 2.25 
(see Liseau et al. 2013). 

References. (1) Hipparcos. (2) Bessell (1990). (3) Engels et al. (1981). 
(4) A. Mora (priv. comm.; FWHM(24yum) = 6"). Binary separation 
on 9 April 2005, 10"4 (5) Hi-GAL: KPOT_smolinar_l, fields 314_0 
& 316_0. Herschel- beams FWHM(10 jjm) = 5'/6, (100 /um) = 6'/8, 
(160 /im) = 11 '.'3, (250 /um) = 17'/6, (350 yum) = 23"9, (500 /im) = 
3572. Binary separation on 21 August 2010, 6'/3. (6) DUNES: 
KPOT_ceiroa_l. Binary separation 29 July 2011, 577. (7) 384. C- 
1025, 380.C-3044(A): FWHM(810/im ) = 19"5. Binary separation 
20-13 November 2007, 8'/8 and 19 September 2009, 770. 

wavelengths beyond 160/rm, the dust emission from the galac- 
tic background is expected to be in the Rayleigh- Jeans (RJ) 
regime, and the emissivity should remain reasonably constant. 
Therefore, when normalised to a specific feature, the background 
is expected to look essentially the same at all wavelengths, when 
disregarding the slight deterioration of the resolution due to the 
smearing at the longest wavelengths. The result of this procedure 
is shown in Fig. 4, where a Cen is represented by the scaled PSFs 
and the background- subtracted fluxes are reported in Table 2. 
Assuming 10-30% higher (lower) fluxes for a Cen results in de- 
pressions (excesses) at the stellar position, which are judged un- 
realistic. This should therefore provide a reasonable estimate of 
the accuracy of these measurements. The photometric calibra- 
tion of SPIRE is described by Bendo et al. (2013). 

2.2. APEX 

The Atacama Pathfinder Experiment (APEX) is a 12 m sub- 
millimetre (submm) telescope located at 5105 m altitude on the 
Llano de Chajnantor in Chile. According to the APEX home 
page 3 , the telescope pointing accuracy is 2" (rms). The gen- 
eral user facilities include four heterodyne receivers within the 
approximate frequency bands 200-1400 GHz and two bolome- 
ter arrays, centred at 345 GHz (870/rni) and 850 GHz (350/mi), 


3 http : //www . apex- telescope . org/telescope/ 


respectively. We used LABOCA and SHeFI APEX-1 for the ob- 
servations of a Centauri. 


2.2.1. LABOCA 870 yum 

The Large APEX BOlometer CAmera (LABOCA) is a submil- 
limetre array mounted on the APEX telescope (Siringo et al. 
2009). The operating wavelength is 870 //m centred on a 150 /i m 
wide window (345 and 60 GHz, respectively). With separations 
of 36", the 295 bolometers yield a circular field of view of 1 1'A. 
The angular resolution is HPBW = 1975 and the under-sampled 
array is filled during the observations with a spiral mapping 
method. 

The mapping observations of a Cen were made during two 
runs, viz. on 10-13 November 2007, and on 19 September 2009. 
The data associated with the programmes 380.C-3044(A) and 
384.C- 1025(A) were retrieved from the ESO archive. Since this 
observing mode is rather inefficient for point sources, we per- 
formed additional test observations with the newly installed tele- 
scope wobbler on 20 May and 13-14 July 2011. While centring 
a Cen B on LABOCA channel 71, the chopping was done by a 
fixed ±25" in east-west direction, which resulted in asymmetric 
sky-flux pickup in this confused field, and these data therefore 
had to be discarded. 

The map data were reduced and calibrated using the software 
package CRUSH 2 developed by AttilaKovacs 4 . The data have 
been smoothed with a Gaussian of HPBW = 13", resulting in an 
effective FWHM = 23 "4. However, fluxes in Jy/beam are given 
for an FWHM = 19"5. 

2.2.2. SHeFI 230 GHz (1300yum) 

Complementing observations were made with the Swedish 
Heterodyne Facility Instrument (SHeFI) APEX-1 during 
Director’s Discretionary Time (DDT) on 16 August 2012. 
APEX-1 is a single- sideband receiver that operates in the 
230 GHz band (21 1-275 GHz) with an IF range of 4-8 GHz. At 
the CO (2-1) frequency, 230.5380 GHz, the half-power-beam- 
width is 27". The estimated main beam efficiency is 7/ m b = 0.75, 
and the kelvin-to-jansky conversion for the antenna is 39 Jy/K at 
this frequency. 

The observing mode was on-the-fly (OTF) mapping with 
a scan speed of 6" s _1 , which resulted in a data cube for the 
5' x5' spectral line map with 441 read-out points (Fig. 6). The 
reference, assumed free of CO emission, was at offset position 
(+1500", -7200"). During the observations, the system noise 
temperature was typically T sys = 200 K. As backend we used 
the fast Fourier transform spectrometer (FFTS) configured to 
8192 channels having velocity resolution 0.16kms _1 , yielding 
a total Doppler bandwidth of 1300 km s _1 . 

These ON-OFF observations generate a data cube with two 
spatial dimensions (the “map”) and one spectral dimension (in- 
tensity vs frequency spectrum). The data were reduced with the 
xs-package by P. Bergman 5 . This included standard, low order, 
base line fitting, and subtraction, yielding for each “pixel” a 
spectrum with the intensity given relative to the zero-kelvin level 
(see Fig. 6a). 


4 See: http://www.submm.caltech.edu/~sharc/crush/ 

download.htm 

5 http : //www . chalmers . se/r ss/oso- en/observations/ 
data-reduction- software 
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Jy/pxl 

0.00113 

0.00091 

0.00068 

0.00046 

0.00024 

0.00002 

- 0.00020 

-0.00043 

-0.00065 




Jy/pxl 

0.0051 

0.0042 

0.0032 

0.0023 

0.0014 

0.00045 

-0.00048 

-0.0014 

-0.0024 


Fig. 2. Binary a Cen at 100 pm (left: green) and 160yum (right: red) at high contrast. Scales are in Jy/pxl, where square pixels are 1" on a side at 
100 pm and 2" on a side at 160 /urn. The field of view is 1 '75 x 3 '5. The image of a Cen is extended in the direction of the companion and, at 100 pm, 
the pair is quasi-resolved, having a separation of 577 along position angle 254°. The white arrows of length 80 AU correspond approximately to 
the stable circumbinary regime around a Cen AB. These point to parts of coherent structures that are also are seen at all FIR/submm wavelengths. 


3. Results 

3.1. Herschel PACS and SPIRE 

The stars were clearly detected at both PACS 100/mi and 
160//m, as expected on the basis of the adopted observing 
philosophy of the DUNES programme. The last entry in Table 1 
refers to the difference between the commanded position, cor- 
rected for proper and orbital motion, and the observed posi- 
tion, determined from the fitting of Gaussians to the data, and 
these are found to clearly be less than 2". This compares very 
favourably with the findings for PACS of the average offset of 
2"4 by Eiroa et al. (2013) for a sample of more than 100 solar 
type stars. Similarly, observed offsets are also acceptable for the 
LAB OCA data, i.e. within 2cr of the claimed pointing accuracy 
(Table 1). 

Observed flux densities of crCenA and B are given in 
Table 2, where the quoted errors are statistical and only refer to 
relative measurement accuracy. Absolute calibration uncertain- 
ties are provided by the respective instrument teams and cited 
here in the text. Together with complementing Spitzer MIPS 
data, these data are displayed in Fig. 7. For both a Cen A and B, 
a marginal excess at 24 pm corresponding to 2.5cr and 2.6cr is 
determined, respectively. The measured flux ratio at this wave- 
length corresponds well with the model ratio, i.e. 2.26 ±0.11 
and 2.25 respectively. 

cr Cen was detected at all SPIRE wavelengths, but with a 
marginal result at 500 pm (Fig. 5 and Table 2). As is also evi- 
dent from these images, the pair was not resolved. 


3.2. APEX LABOCA 

The LABOCA fields are significantly larger than those of the 
PACS frames, and they contain a number of mostly extended 
sources of low intensity. However, cr Cen is clearly detected 
(Fig. 5), but the binary components are not resolved. In Table 2, 
only the average is given for the 2007 and 2009 observations, 
because these flux densities are the same within the errors. 

The proximity of the crCen system leads to angular size 
scales that are rarely ever encountered among debris discs, 
which are generally much farther away. The PACS field of view 
is 105" x 210", so that a Edgeworth- Kuiper belt analogue would 
easily fill the images. The LABOCA frames contain structures 
similar to those discernable in Fig. 2. Although at faint levels, 
most of these sources are definitely real, because they repeatedly 
show up in independent data sets at different wavelengths ob- 
tained with different instruments and at different times. Knotty, 
but seemingly coherent, arcs and ring-like features on arcmin 
scales mimic the morphology of belt features. If these features 
were associated with the a Cen stars, they should move in con- 
cert with them at relatively high speeds. At 870 pm, the strongest 
feature lies less than one arcminute southeast of a Cen and is 
barely detectable at 100 pm, but clearly revealed at 160 pm. 
Measurements of the sky positions for both this feature (“Bright 
Spot”) and trCen are reported in Table 3. For the Bright Spot 
these differ by +172 in right ascension (RA) and +079 in dec- 
lination (Dec.), which provides an estimate of the measure- 
ment error, i.e. ~1" for relatively bright and centrally condensed 
sources. 
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Fig. 3. Simulated observation at 160/rni of crCen AB, assuming a face- 
on circumbinary dust disc/ring {m& ~ 0.1 M ^ , cf. Sect. 4.5.2). As seen 
from the system’s bary centre (white rectangle), the ring extends from 70 
to 105 AU. (The white arrow has a length of 80 AU.) In this image, each 
pixel is 1704 on a side. The image was smoothed with a 16 x 16 pixel 
Gaussian filter with a standard deviation of 3 pixels. The added noise is 
also Gaussian-distributed and smoothed to mimic the noise and beam 
size of Fig. 2. 


Table 3. Dual epoch position data with LABOCA. 


Object 

Year 

RA (J2000) 

Dec (J2000) 

a Centauri AB 

2007 

2009 

14 39 32.174 
14 39 31.246 

-60 50 02.5 
-60 49 59.6 

Bright Spot 

2007 

2009 

14 39 33.953 
14 39 34.120 

-60 50 44.5 
-60 50 43.6 


Owing to its high proper motion, these values are signifi- 
cantly higher for a Cen, viz. A RA = -678 and A Dec = +279, at 
position angle PA ~ 247°. Neglecting the 071 parallactic contri- 
bution, SIMBAD 6 data for a Cen yield for the 1 .856 yr observing 
period A RA = -677, A Dec = +173, and PA = 259°. 

Similar results are obtained involving other features in the 
870 yum images. This provides firm evidence that the inhomoge- 
nous background is stationary and an unlikely part of any cir- 
cumbinary material around crCenAB. A few notes concerning 
the nature of this background are given in the next section. 


6 http : //simbad . u- strasbg . fr/simbad/ 
A102, page 6 of 15 


3.3. APEX-1 SHeFI 

In the direction of trCen, strong CO (2-1) emission lines are 
found at local- standard-of-rest (LSR) velocities of -30, -50, and 
-60kms _1 , with a weaker component also at +40kms _1 . This 
is entirely in accord with the observations of the Milky Way in 
CO (1-0) by Dame et al. (2001). As an example, the distribu- 
tion of the integrated line intensity of the -50kms _1 compo- 
nent is shown in Fig. 6, providing an overview of the molecular 
background in this part of the sky. With the line widths typical 
of giant molecular clouds, this confusing emission is certainly 
galactic in origin and not due to an anonymous IR-galaxy be- 
cause the lines would be too narrow. We can also exclude the 
possibility of a hypothetical circumbinary dust disc around the 
pair a Cen AB, since the observed lines fall at unexpected LSR- 
velocities (v a c en ~ -3kms _1 ) and are also too wide and too 
strong. Since these data do not reveal the continuum, the stars 
a Cen AB are not seen in the CO line maps. 


4. Discussion 

4.1. FIR and submm backgrounds 

At a distance of 1.3 pc and located far below the ecliptic plane, 
any foreground confusion can be safely excluded. At greater 
distances, bright cirrus, dust emission from galactic molecular 
clouds, PDRs, H II regions, and/or background galaxies poten- 
tially contribute to source confusion at FIR and submm wave- 
lengths. However, the projection of crCen close to the galactic 
plane makes background issues rather tricky (see Sect. 2.1.3). In 
particular, extragalactic observers tend to avoid these regions, 
and reliable catalogues for IR galaxies are generally not avail- 
able. However, as our spectral line maps show (Fig. 6), the ob- 
served patchy FIR/submm background is clearly dominated by 
galactic emission. 

4.2. The SEDs of a Cen A and a Cen B 
4.2.1 . The stellar models 

For the interpretation of the observational data we need to ex- 
ploit reliable stellar atmosphere models and stellar physical pa- 
rameters. The model photosphere parameters are based on the 
weighted average as the result of an extensive literature sur- 
vey. The models have been computed by a 3D interpolation 
in a smoothed version of the high-resolution PHOENIX/GAIA 
grid (Brott & Hauschildt 2005) and with the following pa- 
rameters from Torres et al. (2010) and the metallicities from 
Thevenin et al. (2002): (r eff , lo gg, [Fe/H]) = (5824 K, 4.3059, 
+0.195) for a Cen A and (5223 K, 4.5364, +0.231) for orCenB. 
These models are shown in Fig. 7, together with the photom- 
etry (Table 2). We wish to point out again that these data have 
not been used in the analysis, but merely serve to illustrate the 
goodness of the models. 

As a sanity check we compared the integrated model 
fluxes with the definition of the effective temperature, 
fS v dv/fB y (T eff )dv. The luminosities are conserved to within 
0.6% for cr Cen A and 0.1% for crCenB, which is well within 
the observational errors of 2% and 5%, respectively (Torres et al. 
2010), which in turn are comparable to the errors of the theoreti- 
cal model (e.g. Gustafsson et al. 2008; Edvardsson 2008). These 
small differences can probably be traced back to the re-gridding 
of the high-resolution models onto a somewhat sparser spectral 
grid. 
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Fig. 4. Left : at 70 pm with PACS, the binary a Cen AB is clearly resolved (Hi-GAL data). As a result of the re-binning using cubic spline interpo- 
lation, the stars appear too close in this image. Right: the method for the estimation of the source and background fluxes in the SPIRE images is 
illustrated. Shown is a cut in right ascension through the image, normalised to the interstellar dust feature 40" east of a Cen. For comparison, the 
160 pm PACS data are also shown in blue. For SPIRE, green identifies 250 /um, red 350 /um, and black 500 /um. 


4.2.2. The temperature minima of a Cen A and a Cen B 

At heights of some five hundred kilometers above the visible 
photosphere ( h ^ 10 -3 R Q ), the Sun exhibits a temperature min- 
imum, i.e. T /T g ff < 1, beyond which temperatures rise into the 
chromosphere and corona. Theories attempting to explain the 
physics of the heating of these outer atmospheric layers gen- 
erally invoke magnetic fields (Carlsson & Stein 1995), but the 
details are far from understood and constitute an active field of 
solar research (e.g. de la Cruz Rodriguez et al. 2013). Such at- 
mospheric structure can be expected to also be quite common on 
other and, in particular, solar- type stars. There, the dominating 
opacity, viz. H - free-free, limits the visibility to the FIR pho- 
tosphere. For instance, on the Sun, r min occurs at wavelengths 
around 150 /mi. This is very close to the spectral region, where 
this phenomenon has recently been directly measured for the 
first time also on another star (viz. a Cen A, Liseau et al. 2013), 
where in contrast to the Sun, the spatial averaging over the unre- 
solved stellar disc is made directly by the observations. Since it 
is so similar in character, a Cen A may serve as a proxy for the 
Sun as a star (see also Pagano et al. 2004). 

This Tmin -effect has also been statistically observed for a 
larger sample of solar-type stars (see Fig. 6 of Eiroa et al. 2013), 
a fact that could potentially contribute to enhancing the rela- 
tionship between solar physics and the physics of other stars. In 
the end, such observations may help to resolve a long-standing 
puzzle in solar physics concerning the very existence of an ac- 
tual gas temperature inversion at 500 km height and concern- 
ing whether this is occurring in the Quiet Sun or in active re- 
gions on the Sun (e.g., Leenaarts et al. 2011; Beck et al. 2013). 
Theoretical models have hitherto been inconclusive. 

For cr Cen A, Liseau et al. (2013) find that r min = 3920 ± 
375 K, corresponding to the ratio T m i n /r e ff = 0.67 ± 0.06. This 


is lower than observed for the Sun from spectral lines in the 
optical, viz. -0.78 from analysis of the Ca n K-line (Ayres 
et al. 1976). These authors also estimated this ratio for or Cen B, 
i.e. 0.71 and 0.72 using T eff = 5300 and 5150 K, respectively. 
For the data displayed in Fig. 7, we arrive at an estimate of 
Tmin = 3020 ± 850 K near 160 /on for this star. This results in 
Tmin/Teff = 0.58 ± 0.17, with T e ff = 5223 ± 62 K. This value 
for a Cen B is slightly lower than the optically derived value, but 
the error on the infrared ratio is large. However, also for the Sun 
itself, T m i n determinations at different wavelengths show a wide 
spread of several hundred Kelvin (Avrett 2003). 

In the present paper, we are however primarily concerned 
with the possible effects T m in might have on estimation of very 
low emission levels from exo-Edgeworth- Kuiper belt dust. The 
intensity of the stellar model photosphere beyond 20 to 40 pm 
is commonly estimated from the extrapolation of the SED into 
the Rayleigh- Jeans (RJ) regime at the effective temperature T e ff . 
There is a potential risk that this procedure will overestimate 
actual local stellar emissions, which may be suppressed at the 
lower radiation temperatures. In those cases where the SEDs are 
seemingly well fit by the RJ-extrapolations, the differences may 
in fact be due to emission from cold circumstellar dust (exo- 
Edgeworth- Kuiper belts), and here, we wish to quantify the mag- 
nitude of such an effect. 


4.3. Far infrared: detectability of cool dust 
from Edgeworth-Kuiper belt analogues 

The SED of crCenA is exceptionally well-determined over a 
broad range of wavelengths, and we use it here for the estima- 
tion of the likely level of this effect on other, non-resolved sys- 
tems. We do that by “filling the pit with sand” to determine the 
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Fig. 5. Gallery of SPIRE and LABOCA data for a Centauri, which is 
inside the circles. Top to bottom : 250 pm, 350 pm, and 500 pm (SPIRE) 
and 870 yum (LABOCA), where the length of the arrow corresponds 
to 2'. 

corresponding fractional dust luminosity / d = Adust /Tstar and the 
accompanying temperature of the hypothetical dust. 

True Edgeworth- Kuiper belt analogues would be extremely 
faint and very difficult to detect (/ d0 ~ 10 -7 : Teplitz et al. 
1999; Shannon & Wu 2011; Vitense et al. 2012). For estimating 
faint exo-Edgeworth- Kuiper belt levels, we apply simple dust 
ring models to “fill” the flux-dip of the temperature minimum of 
a Cen A, i.e. by minimising the chi-square for 

Xv = , with limiting^ < 10 . (1) 



-100 -50 0 50 
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Fig. 6. a) Grand average of all CO (2-1) spectra, in T mb vs. u L sr, of 
the 5' x 5' map toward a Centauri. Strong emission features at several 
lsr- velocities are evident in this direction of the Galaxy (ft 1 = 315°73, 
b 11 = -0°68), with the strongest line at u LS r = -50kms _1 . b) Integrated 
intensity map of the CO (2-1) line at u LS r = -50kms _1 , i.e. Jr_ 50 dv, 
in the direction of a Centauri (for APEX and at v = 230 GHz, fl dv 
(Kkm s -1 ) = 6 x 10 15 f F v dv (erg cnT 2 s -1 )). The origin of the map is at 
the J2000 equatorial coordinates of a CenB, i.e RA = 14 h 39 m 35 s .06 and 
Dec = -60° 50' 1571, and read-out positions are marked by crosses. 
At the time of observation (16 August 2012), the star is located at the 
centre of the white ring. The bright feature south of the map centre is 
also prominent in the continuum observations long ward of 100 pm (cf. 
Fig. 2) and is the dominating source at 870yum (Fig. 5). 

For illustrative purposes, we place the object at a distance of 
10 pc. Several ring temperatures were tested, and dust fluxes 
were adjusted to match the stellar RJ-SED. The modified black 
body of the dust was then adjusted to coincide with the difference 
in flux (i.e. stellar black body minus measured flux) at 160/mi 
and to a temperature that corresponds to the maximum of this 
flux difference. 

This way of compensating for the T ^ n -effect resulted in a 
dust temperature 7d us t = 53 K at the modified black body radius 
of 34 AU and a formal fractional luminosity /d = (2.2^J 7 ) x 
10 -7 , almost one order of magnitude below the detection limit 
of Herschel- PACS at 100 pm (/ d ~ 10 -6 : see e.g. Fig. 1 of Eiroa 
et al. 2013) and comparable to what is estimated for the solar 
system Edgeworth- Kuiper belt (~10 -7 , Vitense et al. 2012). This 
fit is shown at the bottom of Fig. 8 where it is compared to the 
stellar SED and the difference between the black body and the 
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Fig. 7. SEDs of the binary a Cen AB. The model photospheres for the individual stars are shown by the blue and red lines. The PHOENIX model 
computations extend to about 45 jum (solid lines). At longer wavelengths, these correspond to extrapolations using the proper B v (T e g) (dashes). 
Photometric data points are shown for comparison, with the upper limit being 3<x (Table 2). The green curve shows an ISO-SWS low-resolution 
observation, viz. TDT 60702006 (PI Waelkens; see also Decin et al. 2003). The inset displays the SED of the secondary a Cen B in the FIR/submm 
spectral region. 


temperature minimum fit, whereas the upper panel shows the 
observed data points, the adopted run of the r min - dip and that of 
the stellar RJ-SED. 

Based on simple arguments, this dust temperature also pro- 
vides a first-order estimate of the dust mass emitting at these 
wavelengths (cf. Hildebrand 1983), viz. 


^dust ~ 


[i 


/ m od°lw ] ,s ' model ^ v ^ p2 
^extC 17 ) b(v, r dust ) 


(2) 


For an a Cen-like star not showing the T m i n effect, i.e. where the 
observed S 0 bs,min(v) = S^oddO) (hence M dust = 0), this formula 
would give a “concealed” dust mass of the order of the solar 
Edgeworth-Kuiper Belt (see e.g. Wyatt 2008), i.e. ^10 ~ 3 
for /C 150 //II 1 ^ 10 cm 2 g _1 , depending on the largest size a max , the 
compactness and the degree of ice coating of the grains. For 
the cases considered here, reasonable values are within a fac- 
tor of about 2 to 5 (e.g. Miyake & Nakagawa 1993; Krtigel 
& Siebenmorgen 1994; Ossenkopf & Henning 1994). In gen- 
eral, /CextM = ^scaM + ^absM is the frequency-dependent mass 
extinction coefficient. At long wavelengths, say A > 5 /nm, the 


scattering efficiency becomes negligible for the grains consid- 
ered here, and /c ex t reduces to /c a b S (k v henceforth). 

In the FIR, i.e. for low frequencies, one customarily approx- 
imates the opacity by a power law, i.e. k v = kq (v/vor when 
v < vo, and where To = c/v o > 2 na is a fiducial wavelength in 
the Rayleigh- Jeans regime, e.g. To = 250 yum (Hildebrand 1983). 
For spherical Mie particles, the exponent of the frequency depen- 
dence of the emissivity (3 is in the interval 1 to 2 for most grain 
materials (Emerson 1988). For blackbody radiation, (3 = 0. For 
circumstellar dust, (3 - values around unity have commonly been 
found (e.g. Krtigel & Siebenmorgen 1994; Beckwith et al. 2000; 
Wyatt 2008). 

The rough estimate from Eq. (2) only refers to the a < 
1 mm regime and will thus not account for larger bodies, in 
many cases dominating the mass budget of debris discs (as in 
the solar system). As is commonly done in debris-disc studies, 
we distinguish between the dust mass and the total disc mass. To 
retrieve the mass of the unseen larger objects, numerical col- 
lisional models or analytical laws can be used (e.g. Thebault 
& Augereau 2007; Wyatt 2008; Fohne et al. 2008; Heng & 
Tremaine 20 10a, b; Gaspar et al. 2012), but we focus here solely 
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Fig. 8. Top\ run of the brightness temperature of a Cen A, where the 
solid blue line is a fit through the data points around the temperature 
minimum, and the dashed blue line is the stellar blackbody extrapola- 
tion. The excess beyond 300 pm is presumed to originate in higher chro- 
mospheric layers (e.g. De la Luz et al. 2013). Bottom : SED of a Cen A 
scaled to the distance of 10 pc. The dashed black curve is the difference 
between the black body extrapolation (dashed blue line) and the T min -fit 
(solid blue line). The solid black curve is the emission of the fitted 53 K 
dust ring. 


on the observable dust mass, which can be reasonably derived 
from infrared and subillimetre observations (Heng & Tremaine 
2010a; Heng 2011). Even when obtained with an expression as 
simple as Eq. (2), these estimates are uncertain by at least a fac- 
tor of several, because of uncertainties in the dust opacities. 

From the dynamics of the binary system, we deduce 
(Appendix A) that a circumbinary disc is possible at a distance 
larger than 70-75 AU from the bary centre of the stars (see also 
Wiegert & Holman 1997; Jaime et al. 2012). Comparing the syn- 
thetic image of such a ring with an assumed outer radius of 
105 AU in Fig. 3 with our observed continuum maps (Figs. 2 
and 5), we find seemingly coherent structures at the proper dis- 
tances of a face-on circumbinary ring. It is important to note that 
the circumbinary’ s inclination may be unconstrained in contrast 
to the circumstellar discs that are dynamically limited to incli- 
nations smaller than about 60°, (Wiegert & Holman 1997; see 
also Moutou et al. 2011 for a more general spin-orbit-inclination 
study). Also, Kennedy et al. (2012) report a circumbinary and 
circumpolar dust disc around 99 Herculis, which supports this 
possibility of non-coplanarity. 

Regarding the possible detection of a circumbinary dust disc 
around orCenAB, our single-epoch, multi- wavelength images 
would appear inconclusive. However, both proper motion and 
spectral line data obtained with APEX essentially rule out this 
scenario. If there is really a circumbinary disc or ring around 
a Centauri, it remained undetected by our observations. 


Table 4. Properties of the a Centauri binary. 



a Cen A 

a Cen B 

Sp.type a 
T tfs (K) A 
L st „ (L 0 ) b 
M star ( M 0 A 

Sstar (Re)“ 

Octit (AU) rf 

G2 V 

5824 ± 24 

1 S4Q +0 029 
-0.028 

1.105 ±0.007 
1.224 ±0.003 
2.778 ± 1.476 

K1 V 
5223 ± 62 

0 498+0.025 
u '^^°-0.024 

0.934 ± 0.006 
0.863 ± 0.005 
2.522 ± 1.598 

Common parameters 

Inclination to LOS, i (°) c 

79.20 ± 0.04 

Arg. of periapsis, a> (°) c 

231.65 ±0.08 

Long, of asc. node, £2 (°) c 

204.85 ± 0.08 

Period (yr) 6c 


79.91 + 0.01 

Eccentricity c 


0.5179 + 0.0008 

Distance (pc) 6 


1.348 ±0.035 

Age (yr) e 


(4.85 + 0.50) x 10 9 


References. (a) Kervella et al. (2003); {b) Torres et al. (2010); 
(c) Pourbaix et al. (2002); {d) Holman & Wiegert (1999); (e) Thevenin 
et al. (2002). 

4.4. Mid Infrared: warm zodi-dust in asteroid-like belts 
around a Cen A and a Cen B 

4.4.1. Stable orbits 

The tentative excesses at 24 pm may be due to warm dust. 
However, the binary nature of the aCe n system limits the ex- 
istence of stable orbits to three possibilities. These are one large 
circumbinary disc with a certain inner radius and two circumstel- 
lar discs with maximum (or hereafter, critical) semi-major axes. 
These critical semi-major axes can be found using the semi- 
analytical expression of Holman & Wiegert (1999), viz. 

tfcrit = (c\ + c 2 p + c 3 e + C 4 JU e + c 5 e 2 + c 6 pe 2 } a A b (3) 

where e is the eccentricity of the binary orbit, p = M b /(Ma + 
M b ) i.e. the fractional mass of the stars, «ab is the semi-major 
axis of the binary’s orbit, and the coefficients c\ through ce 
(all with significant error bars) were computed by Holman & 
Wiegert (1999). 

Using known parameters for a Cen (see Table 4), we can de- 
duce that the circumstellar discs can not be larger than 2.78 ± 
1.48 AU around a Cen A and 2.52 + 1.60 AU around crCenB, 
i.e. smaller than ~4 AU around either star (Table 4). 

The putative Earth-mass planet, crCenBb (Dumusque et al. 
2012) is small enough (~1.13 M©) that its Hill radius is just 4 x 
10 _4 AU. Our model discs never reach closer to the star than 

0. 08 AU, and consequently, this planet was neglected in our disc 
modelling. 

Owing to the large errors of the estimated a cr \ u it was use- 
ful to test these limits with simple test-particle simulations of 
mass- and size-less particles. The resulting particle discs will 
also be used as a basis for radiative transfer simulations later. A 
steady state of these are shown in Fig. 9 where the accuracy of 
the Holman & Wiegert (1999) estimates is clear. Each simulation 
run was left for 10 3 orbital periods (i.e. ~8 x 10 4 yr). The dynam- 
ics had earlier been examined also by others (e.g. Benest 1988; 
Wiegert & Holman 1997; Holman & Wiegert 1999; Lis sauer 
et al. 2004; Thebault et al. 2009). 

These size limits are reminiscent of the inner solar system, 

1. e. this opens the possibility of an asteroid belt analogue for 
each star that forms dust discs through the grinding of aster- 
oids and comets. Temperature estimates for the solar system 
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Fig. 9. Face-on circumstellar test-particle discs after ~10 3 periods 
shown with the stars close to periapsis. <2 Cen A is colour-coded blue 
(the left star and its orbit), and a Cen B is colour-coded red (the right 
star and its orbit). The green circles represent a cr i t around the stars and 
the magenta circles show estimates of their respective snow lines. 


zodiacal cloud are around 270 K and that of the fractional lu- 
minosity is about 10 -7 (Fixsen & Dwek 2002; Nesvorny et al. 
2010; Roberge et al. 2012), i.e. /d has a similar magnitude to the 
Edgeworth-Kuiper belt (Vitense et al. 2012). 

4.4.2. Disc temperatures and snow lines 

The Spitzer observations at 24 pm (PSF> 6", Table 2; see 
also Kennedy et al. 2012) would not fully resolve these dy- 
namically allowed discs. Strict lower temperature limits are 
provided by black body radiation from the discs, i.e. T bb = 
V^star/(2R bb ) r eff = 186 K and 147 K, respectively, for a Cen A 
and B, where the ratio of absorption to emission efficiency of 
unity for a black body has been used (e.g. Liseau et al. 2008; 
Lestrade et al. 2012; Heng & Malik 2013). 

The a Cen circumstellar discs are constrained to radii smaller 
than 4AU and should be relatively warm. We need to know, 
therefore, the locations of the snow lines in order to under- 
stand the likelihood of the presence of icy grains on stable orbits 
around these stars. 

The snow line can be defined to be the largest radial dis- 
tance at which the sublimation time scale is larger than all other 
relevant time scales of the system (Artymowicz 1997). Our esti- 
mations are based on the equations and methods by Grigorieva 
et al. (2007) and by Lamy (1974). The sublimation time was 
compared with the orbital time as a lower limit and the sys- 
tem age (4.85 x 10 9 yr) as an upper limit. For the temperature- 
radius relation, we assumed blackbody emitters. This analysis 
resulted in sublimation temperatures between 154 K and 107 K 
for a Cen A, at the radial distances of 4.08 and 8.48 AU, respec- 
tively. Therefore, the snow line of a Cen A seems to be outside 
its dynamically stable region of 2.78 AU, which makes the pres- 
ence of any icy disc grains not very likely (Fig. 9). 

For a Cen B, sublimation temperatures were found to be be- 
tween 157 K and 107 K, which corresponds to 2.23 and 4.81 AU, 
respectively. The lower bound is well inside the critical semi- 
major axis of 2.52 AU (see Fig. 9). Since this is based on an 
initial grain size of 1 mm (smaller sizes move the snow line out- 
wards, while larger ones move it inwards), we can assume that 
the existence of a ring of larger icy grains and planetesimals is 
possible at the outer edge of the a Cen B disc. Such a ring could 


supply the disc with icy grains. However, these cannot be ex- 
pected to survive in such a warm environment for long. Thus 
we assume that the iceless opacities from Miyake & Nakagawa 
(1993) are sufficient in this case. Ossenkopf & Henning (1994) 
have computed opacity models for the growth of ice coatings on 
grains. This conglomeration generally results in higher opacities 
than for their bare initial state. 

We use dust opacities that are based on the models pre- 
sented by Miyake & Nakagawa (1993), but scaled by a gas- 
to-dust mass ratio of one hundred. For a “standard” size dis- 
tribution (Festrade et al. 2012), i.e. n(a) oc a ~ 3,5 , the work by 
Ossenkopf & Henning (1994) gave similar results (see also e.g. 
Krtigel & Siebenmorgen 1994; Stognienko et al. 1995; Beckwith 
et al. 2000). These works had different scopes, addressing specif- 
ically coagulation processes, rather than destructive collisions. 
However, for bare grains these opacites may also be applicable 
to some extent for the dust in debris discs 7 . 

In the following sections we present synthetic SEDs obtained 
with the radiative transfer program RADMC-3D (Dullemond 
2012) to assess the disc configurations in more detail. 

4.5. Energy balance, temperature, and density distributions 

Isothermal blackbody emission from thin rings would not make 
a realistic scenario for the emission from dust belts or discs. In 
particular, we need to specify the run of temperature and den- 
sity, which is obtained from the solution to the energy equation, 
balancing the radiative heating with the cooling by the optically 
thin radiation. 

In this balance equation, the loss term expresses the flux den- 
sity of this thermal dust emission from an ensemble of grains 
received at the Earth, viz. 

Sdust(y)=J' J' 47ra 2 (2v, abs («) n(a,r) da dr 3 (4) 

where 4/ra 2 2v, abs (tf) is the emissivity of the spherical grains with 
radius a, Q v ,^(a) is the absorption coefficient, n(a,r) is the 
volume density of the grains at location r, nB v is the thermal 
emittance at the equilibrium temperature T at r, and D is the 
distance to the source (so that 4 nD 2 f S' dust M dv = Ld us t)- The 
absorption coefficient Q Vy abs is related to the dust opacity through 
Q v abs = (4/3 )K v pa, where p is the grain mass density, <2 v ,abs is 
shown as function of wavelength for a variety of materials and 
compositions for two grain sizes (10 and 1000 pm) in Fig. 4 of 
Krivov et al. (2013). We discuss grain sizes in the next section. 

4.5.1. Grain size distribution 

It is often assumed that n(a) oc a q describes the distribution in 
size of the grains (e.g. Miyake & Nakagawa 1993; Krtigel & 
Siebenmorgen 1994; Ossenkopf & Henning 1994; Krivov et al. 
2000). However, the shape, the existence, and magnitude of the 
power-law exponent q have been a matter of intense debate since 
Dohnanyi’s classic work on interplanetary particles in the so- 
lar system (Dohnanyi 1969). For a collision-dominated system 
at steady state, his result is summarised in a nutshell by the 
parameter v = 11/6, implying that q = 2 - 3s = -3.5 (see 
e.g. Dominik & Decin 2003; Wyatt 2008). This original work 
referred to the zodiacal cloud, but more recent estimates also 
for the much larger “grains”, i.e. Edgeworth-Kuiper belt objects 

7 See also the homepage of B.Draine: http ://www. astro, 
princeton . edu/~draine/dust/dustmix . html 
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(TNOs), by Bernstein et al. (2004) lead to a similar frequency 
spectrum for the diameters, dn/dD oc D~ p and p = 4 + 0.5, so 
that q = -3 ± 0.5 (see also Fraser & Kavelaars 2009). 

More recent works have been following the collisional time 
evolution of dusty debris discs using numerical statistical codes. 
They have found rather strong deviations from true power laws, 
showing wiggles and wavy forms due to resonances (Krivov 
et al. 2006, 2013; Thebault & Augereau 2007; Lohne et al. 2008; 
Krai et al. 2013). However, mean values of these oscillations 
may still be consistent with power-law exponents that are not too 
dissimilar from the steady state Dohnanyi-distribution. In fact, 
Gaspar et al. (2012) find, from a large number of simulations, 
that q ~ -3.65, i.e. slightly steeper than commonly assumed. 
Quoting Miyake & Nakagawa (1993), “q will be large ( — 2.5) if 
the coagulation processes are dominant, whereas q will be small 
( — 3.5) if the disruption processes are dominant”, one might in- 
fer that disruptive collisions would dominate in debris discs, as 
expected. To be consistent with the exploited opacities, but also 
for comparison reasons, we are using a Dohnanyi distribution, 
i.e. n(a) oc a~ 3,5 that we use henceforth. For simplicity, the dis- 
tribution of particle sizes in the disc is assumed to be homoge- 
neous 8 since the grains are large enough to not be significantly 
affected by drag forces and radiation pressure. 

In Eq. (4), the integral over a extends over the finite range 
a m i n to <2 max , where we adopt for a max the upper bound for 
the /^-calculations, i.e. 1 mm. The smallest size, a m \ n , is found 
for grains that are large enough not to be susceptible to the 
combined radiation pressure and stellar wind drag (Strubbe & 
Chiang 2006; Plavchan et al. 2009), i.e. 


^min > ^blow 


8 71 G A/star Pg rain 


+ M st 


( 5 ) 


where we have set the radiative and wind coupling coefficients 
to unity for this order of magnitude estimate. Inside the brackets 
are the radiative and the mechanical momentum rates (“forces”), 
respectively. The stellar mass loss rate of a Cen A is close to that 
of the Sun, i.e. 2 x 10 -14 M© yr -1 , and also the average wind ve- 
locity is similar, viz. i/ w ~ 400 km s -1 (Wood et al. 2001, 2005) 9 . 
This implies that the second term in Eq. (5) is negligible and that 
^biow = 0.64 pm for the stellar parameters of Table 4 and a grain 
density of 2.5 gem -3 (for crCenB, the corresponding a-value is 
about three times lower). It has been shown that the lower cut- 
off size is smooth (e.g. Wyatt et al. 2011; Lohne et al. 2012, and 
references therein) and up to about six times the blow-out radius, 
i.e. a min ~ 4 pm. 

In the next section, we calculate a number of numerical mod- 
els for Eq. (4) by varying parameter values to get a feeling for 
how well (and uniquely) the observations can be reproduced. 
There were marginal excesses corresponding to 2.5<x at 24pm 
for both a Cen A and a Cen B . 


4.5.2. RADMC-3D exozodi-disc calculations 

RADMC-3D is a Monte-Carlo radiative transfer code, written 
primarily for continuum radiation through dusty media in an 

8 In general, grain sizes may not be homogeneously distributed 
throughout circumstellar discs. Smaller grains may be more abundant 
in the outer parts of the discs and may even reside in the dynamical un- 
stable regions in a binary, owing to production through collisions and 
radiation pressure (e.g. Thebault et al. 2010). 

9 In this work the binary is not resolved into its components A and B, 

so that the quoted values refer to both of them. Athough a Cen B is the 

more active star of the two, we simply assumed equal mass loss rates. 


Table 5. Model parameters for RADMC-3D disc calculations. 


Parameter (unit) 

Value 

Lower grain cut-off size, o m i n (pm) 

4 

Upper grain size limit, a max (pm) 

1000 

Power law exponent, n(a) oc a~ q (cm -3 ) 

2.0, 2.5, 3.0, 3.5, 4.0 

Power law exponent, h(r) oc r p (AU) 

0.0, 1.0 

Power law exponent, X(r) oc r~ y (gem -2 ) 

0, 1/2, 1, 3/2, 2, 6, 10 

Grain mass density, p (g cm -3 ) 

2.5 

Vaporization temperature, T va P (K) 

750 

Inner disc radius, 7? in (AU) 

KTvap) 

Outer disc radius, R out (AU) 

^crit 

Logarithmic dust mass (M^ ) 

-7 through +4 by 0.1 


arbitrary three-dimensional geometry (developed by Dullemond 
2012) 10 * * . It might seem an enormous “overkill” to use a radiative 
transfer program for the optically very thin discs discussed here, 
where the optical depth r is comparable to the fractional lumi- 
nosity fa of the dust disc. However, RADMC-3D is a convenient 
and user-friendly means to solve the energy equation, calculate 
SEDs, and convolve the model-images with the appropriate re- 
sponse functions of the used equipment for the observations (e.g. 
transmission curves, PSFs or beams, etc.). Details regarding the 
computational procedure can be found in Appendix A and exam- 
ined parameters are presented in Table 5. 

The basic input to RADMC-3D is a stellar model atmosphere 
(we use 3D interpolations in the PHOENIX/GAIA grid, Brott & 
Hauschildt 2005), a grid with dust densities for each included 
dust species, and absorption and scattering coefficients at each 
wavelength and for each dust species. The inner disc radii are 
given by the dust vaporation temperature 11 of -750 K and the 
outer radii determined by the binary dynamics. The discs have 
assumed height profiles h ~ 0.05 r^, with = 0 (or = 1; e.g., 
Wyatt 2008). 

Also based on the dynamical model parameters are the dust 
mass estimations for various power- law exponents of the surface 
density, X oc r _r , where y > 0. For reference, the solar mini- 
mum mass nebula (SMMN) has y = 1.5 (Hayashi 1981). For the 
SMMN, at r 0 = 1 AU, X 0 ~ 2 x 10 3 gem -2 , whereas exo-planet 
systems may have slightly steeper laws and somewhat higher X 0 
(Kuchner 2004). 

Various grain size distributions were examined in accordance 
with Fig. 7 of Miyake & Nakagawa (1993), i.e. -4 < q < -2. 
The results of these computations are visualised in Fig. 10. 
There, we investigate whether any combination of parameters 
could be consistent with the observed data, i.e. the flux densities 
at 5, 24, and 70//m. For crCenA, acceptable models would im- 
ply values of fa ~ 10 -5 for its hypothetical zodi, where q = -4 
to -3.5. This would appear to agree with the result by Gaspar 
et al. (2012), i.e. q = -3.65. 

However, no such model was found that would fit the data for 
a Cen B, as illustrated in the right-hand panel of Fig. 10. These 
models assume a flat X-distribution, i.e. y = 0. To be compati- 
ble with the observations, models for y ^ 0 would require very 
high values of y (up to 10), making these “discs” more ring-like. 
However, best agreement was still achieved with y = 0, when 
the model disc of a Cen B had less of an extent than the dynam- 
ically allowed size, i.e. Ri n = 0.1 AU and R 0 ut = 0.5 AU. The 
temperature profile is given by T(r) = 745 (r/0.05 AU) -0 55 K, 

10 http://www.ita.uni-heidelberg.de/~dullemond/ 
software/radme- 3d/ 

11 This corresponds approximately to the low-density extrapolation of 
the data for olivine given by Pollack et al. (1994). 
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Fig. 10. Left : mid-infrared SEDs for or Cen A for different values of the power-law exponent q of the size distribution function, viz. n(a ) oc a~ q . 
Black curves depict dust SEDs, whereas blue curves display the total, i.e. S J tar +S J ust . For each q, there are graphs for a selection of initial dust mass 
estimates and the fractional dust luminosity, / d = L dust /L star , is shown next to each curve. The 100 /im data were not part of the fitting procedure. 
Right: same as the left panel but for a CenB and with the total emission shown in red. 



Fig. 11. Mid-infrared SED of a Cen B (star + dust, red curve) with data 
points. Assuming a “standard” disc model with q = 3.5 and y = 1.5, the 
black curve shows the maximum emission from circumstellar dust that 
is consistent with the data. 

and the dust-SED results in an /d = 3 x 10~ 5 for a dust mass 
m d = 4x 10 -6 M(j (3 x 10 20 g, i.e. about 10-30 m Zo ) for grains 
with a = 4//m to 1 mm and with sizes according to a q = -3.5 
power law (see Fig. 11). 


For both a Cen A and B, these limiting values of fd and 
respectively, would correspond to some 10 2 and 10 times those 
of the solar system zodi dust in ^100 //m size particles (Roberge 
et al. 2012; Fixsen & Dwek 2002; Nesvorny et al. 2010). Fixsen 
& Dwek (2002) estimate the total mass of the asteroid belt at 
about 0.03 M(£ (2 x 10 24 g). Heng & Tremaine (20 10a, b) have 
calculated total disc masses for ages of about 3 Gyr. Their Fig. 7 
would indicate that the mass of a “collision-limited” disc with 
a semi-major axis ^4 AU and detected at 24 jum would amount 
to ~0.1 M(£ . In both cases, these estimations refer effectively 
to discs around single stars, and it is not entirely clear to what 
degree this also applies to multiple systems. 

5. Conclusions 

Based on FIR and submm photometric imaging observations 
from space (Herschel) and the ground (APEX) of a Cen, the fol- 
lowing main conclusions can briefly be summarised as 

• Both binary components A and B were detected individually 
at 70/rm. At longer wavelengths, 100-870/mi, the pair was 
unresolved spatially during the time of observation, i.e. 8'.'4 
in 2007 to 5"7 in 2011. 
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• At 100 and 160yum, the solar twin arCenA displays a tem- 
perature minimum in its SED, where the radiation originates 
at the bottom of its chromosphere. At longer wavelengths 
higher radiation temperatures are observed. This T min -effect 
is also marginally observed for crCenB. 

• For solar-type stars not displaying the T m i n -phenomenon, 
compensating for the “missing flux” by radiation from op- 
tically thin dust could account for fractional dust luminosi- 
ties, /d ~ 2 x 10~ 7 , comparable to that of the solar system 
Edgeworth- Kuiper belt; 

• Combining Herschel data with mid-infrared observations 
C Spitzer ) indicates marginal (2.5cr) excess emission at 24 fim 
in the SEDs of both stars. If caused by circumstellar emission 
from dust discs, fractional luminosity and dust mass levels 
would be some 10 to 100 times those of the solar zodiacal 
cloud. 
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star; i.e., they had 10 4 particles in discs with constant particle 
density and radii smaller than 5 AU. The circumbinary disc also 
had constant particle density but was 100 AU in radius, (This 
was varied, but we only show results from the disc with 100 AU 
radius.) The circumbinary disc’s outer radius is, in contrast to the 
circumstellar discs, not constrained by the employed dynamics. 
Therefore the outer edge seen in Fig. 3 is artificial. 

All discs were left to run for 10 3 periods, i.e. ~8 x 10 4 yr, 
which was found to be sufficient to achieve a steady state. That 
the disc had settled by that time was verified with much longer 
runs (e.g. ~8 x 10 5 yr). At the inner edge of the circumstel- 
lar discs, the grid resolution was up to 0.047 AU. On the much 
larger scale of the circumbinary disc, time steps of 0.1 yr were 
sufficient. 

Finally, particles that came either closer than 0.001 Hill radii 
to each star or farther out than 1000 AU from the bary centre were 
removed. These limits were arbitrarily set on the basis of previ- 
ous tests. The number of particles remaining varied from run to 
run with e.g. 2496 as the lowest (crCenB with initial disc ra- 
dius of 5AU) up to 6717 (again orCenB but with initial disc 
radius of 0.2 AU). The crCen A disc with initial radius of 5 AU 
had 3538 remaining particles. 

To compute dust SEDs, the dust disc density grid was ob- 
tained directly from the steady state discs. These were con- 
verted to a grid of mass densities with a large number of initially 
guessed dust masses, ra^st* The proper temperature and density 
profiles yielded SEDs, which were compared to the data. The 
best fit was selected on the basis of its x 1 - value. 


Appendix A: Orbital dynamics and disc modelling References 


The test particle simulations were done for two reasons: firstly, 
to assess the binary dynamics per se and, secondly, for the steady 
state results to be used as disc models for radiative transfer 
simulations. 

The standard equation of motion for circumstellar dust grains 
(around a single star) sums up most of the physical processes a 
grain experiences and is written as (Robertson 1937; Stark & 
Kuchner 2008) 


GM* 


(1 -S„)r„ 


( l+SW)S n GMi, 


( r n r + v n ) 


V* G mi 


(A.I) 


for each test particle n. Here & n is the ratio between radiation 
pressure and gravitation exhibited on particle n, S W is the ratio 
between stellar-wind drag and Poynting-Robertson drag and is 
for the Sun about one third (~0.2-0.3, Gustafson 1994), and m* 
and r t are the mass and position of any planet in the system. 

As hinted in Sect. 4.5.1, we consider grain sizes >4 /mi. This 
is more than six times the blow-out radius, and consequently, we 
only solve the classical restricted three-body-problem, i.e. 


a n 


GM a 

Vn ~ '•aI 3 


(r„ - r A ) - 


GM _ b 

Vn - J"b | 3 


(r„ - r B ) 


(A. 2) 


for each test particle n. These simulations implemented a 
Runge-Kutta 4 integrator to solve the motions. 

One particle disc was simulated at a time. The two circum- 
stellar discs had identical initial conditions except for the central 
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